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Abstract 

The coupled cluster method is applied to a spin-half model at zero temperature (T = 0), which 
interpolates between Heisenberg antiferromagnets (HAF's) on a kagome and a square lattice. With 
respect to an underlying triangular lattice the strengths of the Heisenberg bonds joining the nearest- 
neighbor (NN) kagome sites are Ji > along two of the equivalent directions and J2 > along the 
third. Sites connected by J2 bonds are themselves connected to the missing NN non-kagome sites 
of the triangular lattice by bonds of strength J[ > 0. When J[ = Ji and J2 = the model reduces 
to the square-lattice HAF. The magnetic ordering of the system is investigated and its T = phase 
diagram discussed. Results for the kagome HAF limit are among the best available. 

PACS numbers: 75.10.Jm, 75.30.Gw, 75.40.-s, 75.50.Ee 



I. INTRODUCTION 

Quantum magnets defined on two-dimensional (2D) spin lattices exhibit a wide range of 
physical states at zero temperature, from those with classical-type ordering (albeit reduced 
by quantum fluctuations) to valence-bond solids and spin liquids. ^'^ The behavior of these 
strongly correlated and often highly frustrated systems is driven by the nature of the un- 
derlying crystallographic lattice, by the number and range of the magnetic bonds, and by 
the spin quantum numbers of the atoms localized to the lattice sites. Very few exact results 
exist for such 2D systems and the application of approximate methods has become crucial to 
their understanding. A complete picture of their behavior has only slowly begun to emerge 
by considering a wide range of possible scenarios in related models that are themselves often 
inspired, or followed closely afterwards, by their experimental realisation and study. For 
easy and accurate comparisons to be made it is clearly preferable to use the same theo- 
retical technique. Among the most accurate, most universally applicable, and most widely 
applied to quantum magnets of such methods is the coupled cluster method (CCM).^"^ Our 
aim here is to use the COM to extend our understanding of frustrated quantum magnets by 
applying it to a novel 2D system that, in some well-defined sense described below, interpo- 
lates between Heisenberg antiferromagnets (HAF's) defined on square and kagome lattices 
respectively. 

An archetypal and much studied model in quantum magnetism is the frustrated spin- 
half J1-J2 HAF model on the square lattice with nearest-neighbor (NN) bonds (of strength 
Ji > 0) competing with next-nearest-neighbor (NNN) bonds (of strength J2 = aJi > 0). It 
exhibits two different quasiclassical phases with collinear magnetic long-range order (LRO) 
at small {a < Uc^ ~ 0.4) and large {a > ac^ ~ 0.6) values of a, separated by an intermediate 
quantum paramagnetic phase with no magnetic LRO in the regime ttci < a < otc^- Interest 
in this model has been reinvigorated of late by its experimental realisation in such layered 
magnetic materials as Li2VOSi04,^'^ Li2VOGe04,^ VOMo04,^ and BaCdVO(P04)2.^ The 
syntheses of such layered quasi-2D materials has stimulated a great deal of renewed interest 
in the model (and see, e.g., (Refs. [10-13]). Amongst several methods applied to the J1-J2 
model has been the CCM.i^-^^ 

In view of the huge interest in the J1-J2 model there have been several recent attempts to 
investigate various generalizations and modifications of the model, in order to shed further 



light on its properties. As an example of a generalization of the model we mention a recent 
study^^ of the effects of interlayer couplings on the 2D J1-J2 model. This study also employed 
the CCM in its analysis. Modifications of the J1-J2 model that have been studied include 
models wherein some of the NNN J2 bonds are removed. Various such models exist in which 
either half or three-quarters of the J2 bonds are removed in particular arrangements, as 
discussed below. All of these models studied to date have fascinating magnetic properties 
and ground-state phases in their own right. 

One such model is the spin-half anisotropic HAF on the 2D triangular lattice, which 
has also been studied by the CCM,^° and which interpolates between HAF's on square and 
triangular lattices. It is fully equivalent to a variant of the square-lattice J1-J2 model in 
which half of the J2 bonds are removed, leaving just one NNN bond across the same diagonal 
of each basic square plaquette. Thus, for this model, the two cases J2 = and J2 ^> 00 
relate to a HAF on the square lattice and a set of decoupled one-dimensional HAF chains 
respectively, with the HAF on the triangular lattice in between at J2 = Ji- Strong evidence 
was found^° that quantum fiuctuations for this spin-half model favor a weakly first-order (or 
possibly second-order) transition from Neel order to a helical state at a first critical point 
at cici = 0.80 ±0.01 by contrast with the corresponding second-order transition between the 
equivalent classical states at aci = 0.5. The CCM was also, uniquely, powerful enough to 
provide strong evidence for a second quantum critical point at ac2 = 1-8 ± 0.4 where a first- 
order transition occurs between the helical phase and a collinear stripe-ordered phase with 
no classical counterpart, thereby providing quantitative verification of an earlier qualitative 
prediction of such a transition from a renormalization group analysis of the model. ^^ (As a 
parenthetical note, a different way of removing half the J2 bonds, results in the so-called 
Union Jack model, to which the CCM has also recently been applied, ^^ and which has a 
quite different zero-temperatue phase diagram). 

A further modification of the original spin-half J1-J2 niodel is now to remove another 
half of the J2 bonds, leaving half the fundamental square plaquettes with one J2 bond 
and the other half with none. One way of doing this in a regular fashion results in the 
Shastry-Sutherland modeP^ in which no J2 bonds meet at any lattice site and every site 
is five-connected (by four NN Ji bonds and one J2 bond). Interest in this model has been 
renewed by the discovery of the magnetic material SrCu2 (603)2 [Ref. 24] that can be 
understood in terms of it. Its classical ground state is the collinear Neel state for J2/J1 < 1 



and a noncollinear spiral for J2/ Ji > 1, with a second-order phase transition in between. 
However, it is known that the spin-half model has a quantum ground state which is a product 
of local pair singlets (the so-called orthogonal-dimer state) for J2/J1 > 1.465 ±0.025, which 
has no classical counterpart. The CCM has also been applied to this modeP^'^^ and the 
latest results^^ strongly suggest that no intermediate phase exists between the Neel and 
dimerized phases, and that the direct transition between them is of first-order type. 

A different, but equally important, archetypal magnetic system showing frustration, but 
now of the geometric kind rather than the dynamic kind, is the spin-half kagome-lattice HAF. 
Although this system has been the subject of intense study over a long period, the nature 
of its gound state is still not definitively settled. Among the leading theoretical contenders 
are a valence-bond solid state^^"^^ and a spin-liquid state. ^^"^'^ The spin-half kagome-lattice 
HAF has become the subject of renewed interest after a possible physical realization of the 
model has been found experimentally in the herbertsmithite material ZnCu3(OH)6Cl2.^^'^^ 
A spatially anisotropic version of the spin-half kagome-lattice HAF has also been experimen- 
tally studied after its physical realization in the volborthite material Cu3V207(OH)2-2H20.^^ 
This latter model has also been studied theoretically in recent years. ^^"^"^ 

In the present paper we investigate the phase diagram of a novel spin-half HAF that 
is another depleted modification of the J1-J2 model, and which also contains both the 
spatially anisotropic and the isotropic kagome-lattice HAF's discussed above as limiting 
cases. As described in more detail below the model also interpolates continuously between 
the geometrically frustrated kagome-lattice HAF and the unfrustrated square-lattice HAF. 
After describing the model in Sec. II, we apply the CCM to investigate its ground-state 
properties. The CCM is first described briefiy in Sec. HI, and the results are presented in 
Sec. IV. We conclude in Sec. V with a discussion of the results and a comparison of them 
with other results for limiting cases of our model. 

II. THE MODEL 

In this paper we consider an alternate and novel variant of the J1-J2 model in which 
three-quarters of the J2 bonds are removed from the original J1-J2 niodel as in the Shastry- 
Sutherland model above, but in a different pattern, as shown in Fig. 1. It may equivalently 
be obtained from the spin-half anisotropic HAF on the 2D triangular lattice by removing 
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FIG. 1: (Color online) The interpolating kagome-square model; — Ji; - - J(; - ■ - J2; 
showing (a) the canted state; and (b) the semi-striped state. 



every second line of NNN J2 bonds. The square-lattice representation of the model shown 
in Fig. 1 contains the two square sublattices of A sites and B sites respectively, and each 
of these in turn contains the two square sublattices of Ai and A2 sites, and Bi and B2 
sites respectively, as shown. It is very illuminating to consider the anisotropic variant in 
which half of the Ji bonds are allowed to have the strengths J[ > along alternating rows 
and columns. All of the bonds joining sites i and j are of standard Heisenberg type, i.e., 
proportional to Sj-Sj, where the operators Sj = (sf,s^,sf) are the quantum spin operators 
on lattice site i, with s^ = s{s + 1) and s = | for the quantum case considered here. 

This so-called interpolating kagome-square model differs principally from the Shastry- 
Sunderland model in that the A sites are six-connected (by 2 NN Ji bonds, 2 NN J[ bonds 
and 2 NNN J2 bonds) while the B sites are four-connected (by four NN J[ bonds for the 
Bi sites and four NN Ji bonds for the B2 sites). The spin-half HAF's on the 2D kagome 
and square lattices are represented respectively by the limiting cases { Ji = J2, Ji = 0} 
and { Ji = J[, J2 = 0}. The limiting case when {Ji = J[ = 0; J2 > 0} represents a set of 
uncoupled ID HAF chains. The case J[ = with J2 7^ Ji represents a spatially anisotropic 
kagome HAF considered recently by other authors, ^^"^° especially in the quasi-lD limit 
where J2/J1 ^ 1-^° Henceforth we set Ji = 1 and consider the case when all bonds are 
antiferromagnetic, (i.e., J[ > 0, J2 > 0). 

Considered as a classical model (corresponding to the case where the spin quantum num- 
ber s — >• 00) the interpolating kagome-square model has only two ground-state (gs) phases 



separated by a continuous (second-order) phase transition at J2 = J 2 = \{.Ji + J'l)- For 
J2 < Jf the system is Neel-ordered on the square lattice, while for J2 > J2 the system 
has noncollinear canted order as shown in Fig. 1(a), in which the spins on each of the Ai 
and the A2 sites are canted respectively at angles (tt =1= 0) with respect to those on the B 
sublattice, all of the latter of which point in the same direction. The lowest-energy state in 
the canted phase is obtained with = 0ci = cos^"'^(J|'/J2). The Neel state, for J2 < J2, 
simply corresponds to the case 0ci = 0. When J[ = and J2 = Ji, corresponding to the 
isotropic kagome-lattice HAF, 0ci = l^r, as required by symmetry. We also note that as 
J2 ^ 00 (with Ji and J[ finite), 0ci -^ \'^-, and the spins on the A sublattice become antifer- 
romagnetically ordered, as is expected, and these spins are orientated at 90° to those on the 
ferromagnetically-ordered B sublattice. Of course there is complete degeneracy at this clas- 
sical level in this limit between all states for which the relative ordering directions for spins 
on the A and B sublattices are arbitrary. The spin-half problem in the same limit should 
also comprise decoupled ferromagnetic and antiferromagnetic sublattices. We expect that 
this degeneracy in relative orientation might be lifted by quantum fiuctuations by the well- 
known phenomenon of order-hy- disorder. ^"^ Since it is also true that quantum fiuctuations 
generally favor collinear ordering, a preferred state is thus likely to be the so-called ferri- 
magnetic semi-striped state shown in Fig. 1(b) where the A sublattice is now Neel-ordered 
in the same direction as the B sublattice is ferromagnetically ordered. Alternate rows (and 
columns) are thus ferromagnetically and antiferromagnetically ordered in the same direction 
in the semi-striped state. 

III. COUPLED CLUSTER METHOD 

The COM (see, e.g., Refs. [3-5] and references cited therein) that we employ here is one of 
the most powerful and most versatile modern techniques available to us in quantum many- 
body theory. It has been applied very successfully to various quantum magnets (see Refs. [4, 
5,14-20,25,26,52,53] and references cited therein). The method is particularly appropriate for 
studying frustrated systems, for which some of the main alternative methods either cannot 
be applied or are sometimes only of limited usefulness, as explained below. For example, 
quantum Monte Carlo (QMC) techniques are particularly plagued by the sign problem 
for such systems, and the exact diagnoalization (ED) method is restricted in practice by 



available computational power, particularly for s > 1/2, to such small lattices that it is 
often insensitive to the details of any subtle phase order present. 

The method of applying the CCM to quantum magnets has been described in detail 
elsewhere (see, e.g., Refs. [4,5,14,19,52,53] and references cited therein). It relies on building 
multispin correlations on top of a chosen gs model state |$) in a systematic hierarchy 
of LSUBn approximations (described below) for the correlation operators S and S that 
parametrize the exact gs ket and bra wave functions of the system respectively as |^) = e^\^) 
and (^1 = (<l>|S'e^'^. In the work presented here we use two different choices for the model 
state 1$), namely the classical antiferromagnetic Neel state and the ferrimagnetic canted 
state. We note that the ferrimagnetic semi-striped state provides another possible choice 
of model state 1$), but we do not consider it further in the present paper, except in brief 
remarks at the end of Sec. IV. 

In each case we employ the well-established LSUBn approximation scheme in which all 
possible multi-spin-flip correlations over different locales on the (square) lattice defined by n 
or few contiguous lattice sites are retained. As usual the number of independent fundamental 
clusters (i.e., those that are inequivalent under the symmetries of the Hamiltonian and of 
the model state) increases rapidly with the truncation index n. For example, the number 
of such fundamental clusters for the canted model state is 201481 at the LSUB8 level of 
approximation in the triangular-lattice geomery where J2 bonds are considered to join NN 
pairs, and this is the highest level for the present model that we have been able to attain with 
available computing power. In order to solve the corresponding coupled sets of CCM bra- 
and ket-state equations we use massively parallel computing,^^ typically using 600 processors 
simultaneously. We present results below both at various LSUBn levels of approximation 
with n = {2, 4, 6, 8} for the case J[ = Ji = 1 and with n = {2, 4, 6} for other values of the 
bond strengths, and at the corresponding n — > 00 extrapolation (LSUBoo) based on the 
well-tested extrapolation schemes described below and in more detail elsewhere. "^'^''^'^"^^ We 
note that, as always, the CCM exactly obeys the Goldstone linked-cluster theorem at every 
LSUBn level of approximation. Hence we work from the outset in the limit A^ ^ 00, where 
N is the number of sites on the square lattice, and extensive quantities like the gs energy 
are hence always guaranteed to be linearly proportional to A^ in this limit. We note for later 
purposes that the number of sites on the kagome lattice (obtained by removing all Bi sites 
in Fig. 1) is clearly Nx = jN. 



TABLE I: Number of fundamental LSUBn configurations (Nf) for the semi-striped and 
canted states of the spin-1/2 Ji-J[-J2 interpolating square-kagome model. 
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LSUB4 


32 


200 


LSUB6 


645 


6041 


LSUB8 


14936 


201481 



Note that for the canted phase we perform calculations for arbitrary canting angle 
shown in Fig. 1(a), and then minimize the corresponding LSUBn approximation for the en- 
ergy -ELSUBn(0) with respect to to yield the corresponding approximation to the quantum 
canting angle 0LSUBn- Generally (for n > 2) the minimization must be carried out computa- 
tionally in an iterative procedure, and for the highest values of n that we use here the use of 
supercomputing resources was essential. Results for the canting angle 0LSUBn will be given 
later. 

As always, we choose local spin coordinates on each site for each choice of model state, so 
that all spins in |$), whatever the choice, point in the negative z-direction (i.e., downwards) 
by definition in these local coordinates. Then, in the LSUBn approximation all possible 
multi-spin-flip correlations over different locales on the lattice defined by n or fewer contigu- 
ous lattice sites are retained. The operator S thus contains only linear sums of products 
of creation operator s^ = s^ + is\ on various sites A;, while the operator S contains only 
similar linear sums of products of destruction operators s~j^ ^ s^ — is\. The numbers Nf 
of such distinct (i.e., under the symmetries of the lattice and the model state) fundamental 
configurations of the current model in various LSUBn approximations are shown in Table L 
We note that the distinct configurations given in Table I are defined with respect to the 
geometry described in Sec. II, and in which the B sublattice sites of Fig. 1(a) are defined 
to have four NN sites joined to them by either Ji bonds or J[ bonds, and the A sublattice 
sites are defined to have the six NN sites joined to them by Ji, J(, or J2 bonds. If we had 
chosen instead to work in the square-lattice geometry every site would have four NN sites. 

A significant extra computational burden arises here for the canted state due to the 



need to optimize the quantum canting angle at each LSUBn level of approximation, as 
described above. Furthermore, for many model states the quantum number s'^ = X]i=i -^f ^^ 
the original global spin-coordinate frame, may be used to restrict the numbers of fundamental 
multi-spin-flip configurations to those clusters that preserve sf. as a good quantum number. 
This is true for the Neel state where s^ = and for the semi-striped state for which 
s^ = ^Ns, where A^ is the number of lattice sites. However, for the canted model state that 
symmetry is absent, which largely explains the significantly greater number of fundamental 
configurations shown in Table I for the canted state at a given LSUBn order. Hence, the 
maximum LSUBn level that we can reach here for the canted state, even with massive 
parallelization and the use of supercomputing resources, is LSUB8. For example, to obtain 
a single data point for a given value of J2, with Ji = 1 and J[ = 1, for the canted phase 
at the LSUB8 level typically required about 0.2-2.0 h computing time using 600 processors 
simultaneously for non-critical regions. However, for values of J2 near to critical points, 
the LSUB8 computing time increased significantly, typically to lie in the range of 5-24 h to 
obtain a single data point using 600 processors simultaneously. 

At each level of approximation we may then calculate a corresponding estimate of the gs 
expectation value of any physical observable such as the energy E and the magnetic order 
parameter, M = — -^ X]i=i(^kll^)) defined in the local, rotated spin axes, and which thus 
represents the average on-site magnetization. Note that M is just the usual sublattice (or 
staggered) magnetization per site for the case of the Neel state as the CCM model state, for 
example. 

It is important to note that we never need to perform any finite-size scaling, since all 
CCM approximations are automatically performed from the outset in the infinite-lattice 
limit, N -^ 00, where A^ is the number of lattice sites. However, we do need as a last 
step to extrapolate to the exact n ^ 00 limit in the LSUBn truncation index n, at which 
the complete (infinite) Hilbert space is reached. We use here the well-tested^^'^^ empirical 
scaling laws 

E/N = ttQ + airr'^ + a2n~* , (1) 

M = 60 + hn'^ + 62^"^ . (2) 
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FIG. 2: (Color online) Ground-state energy per spin of the spin-1/2 interpolating 

kagome- square model with Ji = J( = 1, using the LSUB4 approximation of the CCM with 

the canted model state, versus the canting angle 0. For J2 < 1.392 the minimum is at 

= (Neel order) at this level of approximation, whereas for J2 ^ 1.392 the minimum 

occurs at = 0lsub4 7^ 0, indicating a phase transition at J2 ~ 1.392 in this 

approximation. Results are shown for those values of for which the corresponding CCM 

equations have real solutions. 

IV. RESULTS 

We show in Fig. 2 our CCM results, based on the canted model state, for the gs energy 
per spin, E/N, plotted as a function of the canting angle 0. We show results specifically at 
the LSUB4 level of approximation for the case Ji = J( = 1, but our results are qualitatively 
similar at other LSUBra levels and for other values of J[ (with Ji = 1). Curves like those 
in Fig. 2 show that at this LSUB4 level of approximation, with Ji = J[ = 1, the minimum 
energy is at = for J2 < J^'^^^'^ ^ 1.392 and at a value 7^ for J2 > J^'^^^'^. Thus, we 
have a clear indication of a shift of the critical point at J2 = J2^ between the quantum Neel 
and canted phases from the classical value J^' = 1 when Ji = J[ = 1. The observation that 
Neel order survives beyond the classically stable regime, for the quantum spin-half system, 
is an example of the promotion of collinear order by quantum fluctuations, a phenomenon 
that has been observed in many other systems. 

In Fig. 3 we show the canting angle 0LsuBn that minimizes the gs energy -ELSUBn(0) at 
various CCM LSUBn levels based on the canted state as model state, with n = {2,4,6,8}, 
again for the case J[ = Ji = 1. We see clearly that at each LSUBn level shown there is a 
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FIG. 3: (Color online) The angle 0LSUBn that minimizes the energy -ELSUBn(0) of the 

spin-half interpolating kagome-square model with Ji = J[ = 1, versus J2- The LSUBn 

approximations with n = {2, 4, 6, 8}, using the canted model state, are shown. The 

corresponding classical result (pci is shown for comparison. 

finite jump in 0LSUBn at the corresponding LSUBn approximation for the phase transition 
at J2 = J^"^^^" between the Neel state (with 0LSUBn = 0) and the canted state (with 
0LSUBn 7^ 0). Thus at each LSUBn level of approximation the quantum phase transition 
is first-order, compared to the second-order classical counterpart. A close inspection of the 
insert in Fig. 3 shows that we cannot completely rule out as n — i> 00, with increasing level 
of LSUBn approximation, the possibility that the phase transition at J2 = J; 
becomes of second-order type, although a weakly first-order one seems more likely on the 
evidence so far. We note that the LSUBn estimates for the phase transition between Neel 
and canted phases, J^^^^", fit well to an extrapolation scheme J^^^^*^ = jlsuboo _|_ ^^-i_ 
The corresponding estimates for Jg^ = jlsuboo ^^ ^j^g ^^^^ shown (Ji = J[ = 1) are Jg^ = 
L298 ± 0.003 based on n = {2, 4, 6, 8} and J^' = 1.302 ± 0.001 based on n = {2, 4, 6} where 
the errors quoted are those associated with the specific fit determined via least-squares. We 
also see from Fig. 3 that results for the quantum canting angle converge very rapidly with 
increasing LSUBn level of approximation, except for a small region near the phase transition 
at J2 = J2^. We note too that as J2 -^ 00 the canting angle (p ^ ^n considerably faster 
than does the classical analog 0ci- Similar estimates for Jg^ have been calculated for other 
values of J[ (with Ji = 1). Thus in Fig. 4 we compare the phase boundary between the 
Neel and canted states for the present spin-half model with its classical analog. We note 
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FIG. 4: Ground-state phase diagram of the spin-half interpolating square-kagonie lattice 
model, calculated as described in the text, compared with its classical counterpart. 
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FIG. 5: (Color online) Ground-state energy per spin versus J2 for the Neel and canted 
phases of the spin-1/2 interpolating kagome-square model (Ji = 1). The CCM LSUBoo 

results based on LSUBn approximations with n = {2, 4, 6, 8} for J[ = 1.0 and with 

n = {2,4,6} for the other two J[ values, are shown. The CCM results using the canted 

model state are based on extrapolated LSUBn approximations with the canting angle 

(p = 0LSUBn that minimizes Ei^suBn{4>)- 



in particular that at J( = the critical value calculated as above is at J2^ = 0.51 ± 0.01, 
compared with the classical value J2 = 0.5 at J[ = 0. 

Our CCM results for the gs energy per spin, E/N, are shown in Fig. 5 as a function 
of J2, for various values of J[ (with Ji = 1). (Note that for the case J[ = 0, we typically 
use a very small value J[ ~ 10^^.) At the isotropic kagome point (J( = 0, J2 = Ji = 
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FIG. 6: (Color online) Ground-state magnetic order parameter versus J2 for the Neel and 

canted phases of the spin-half interpolating kagome-square model (Ji = 1). The CCM 

LSUBcx) results based on LSUBn approximations with n = {2, 4, 6, 8} for J[ = 1.0 and 

■with n = {2, 4, 6} for the other J[ values, are shown. The CCM results using the canted 

model state are based on extrapolated LSUBn approximations with the canting angle 

= 0LSUBn that miniminizes -ELSUBn(0)- 

1) the gs energy per spin is E/N ^ —0.324 ± 0.002, where the error estimate is based 
on comparing LSUBn extrapolations from the three sets n = {2,4,6,8}, {2,4,6}, and 
{4,6,8}. Expressed equivalently in terms of the number N^ of kagome sites, this result 
is E/Nk ~ —0.432 ± 0.002. Our corresponding result for the pure square-lattice HAF 
(with J[ = Ji = 1, J2 = 0) is E/N ^ -0.6697 ± 0.0003. One observes weak signals of a 
discontinuity in the first derivative of the energy at the phase transition points 32^ for all 
values of J(. 

Much clearer evidence of a first-order phase transition at J2^ is seen in Fig. 6 where 
we show comparable CCM results for the average on-site magnetization or magnetic or- 
der parameter, M defined in Sec. III. From the symmetry of the model under the 
interchange of the bonds Ji ^ J(, we note that the order parameter should satisfy 
the relation M{Ji, J[, J2) = M(J(,Ji,J2). Since the order parameter is independent 
of an overall scaling of the Hamiltonian, we may also express this relation in the form 
M(l, J[/Ji, J2/J1) = M{1, Ji/J[, J2/ J'l)- Thus, placing Ji = 1, as we have done here, and 
considering M = M{J[; J2) as a function of the remaining two parameters, we have the 
exact relation M{J[] J2) = M(l/J(; J2/ Jij- The curves shown in Fig. 6(a) for J[ = 0.25 
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and 4.0 and for 0.5 and 2.0 are easily seen to satisfy this relation for both the Neel and 
canted phases, thereby providing a very good check on our numerics. 

For the anisotropic kagome lattice {J[ = 0), the minimum LSUBoo value of M ^ 0.145 
is seen from Fig. 6(b) to occur precisely at the isotropic kagome HAF point {J[ = 0, J2 = 
Ji = 1). We may easily re-express this at the kagome point J( = in terms of M^ = 
— j^ '^i=i{Si), where the sum is taken only over the kagome-lattice sites and where again 
the spins are defined in the local, rotated spin axes in which all spins in the CCM model state 
point in the negative z-direction. Thus, the non-kagome spins on the Bi sites are then frozen 
in the case J( = to have their spins exactly aligned along the local z-axes, and hence at the 
kagome point M^ = ^{M — 0.125) ~ 0.026. Our result is thus that only 5% of the classical 
ordering remains for the spin-half kagome HAF, with an error that makes this compatible 
with zero. Our corresponding result for the square-lattice HAF is M = 0.310 ± 0.003. 

For the Neel phase (0 = 0) curves in Fig. 6 we show the extrapolated (LSUBoo) results for 
all values of J2 (and given values of Ji and J[) for which M > 0, and hence these extend into 
a regime that is unphysical for the Neel phase since the canted phase has lower energy there. 
By contrast, for the canted phase (0 7^ 0) curves in Fig. 6 we show only the extrapolated 
(LSUBoo) results for regimes of J2 (and given values of Ji and J[) for which we have LSUBn 
data (with 0LSUBn 7^ 0) for all of the set n = {2, 4, 6}, in order that the n ^ 00 extrapolation 
can be robustly performed. Curves like those in Fig. 3 thus show that the results for the 
canted phase shown in Fig. 6 terminate (artificially) at J^^^^^ > ■J2 °° = J2^, rather 
than at the physical value J2^. For this (unphysical, but computationally imposed) reason, 
corresponding pairs of curves for M versus J2 (for the same given values of Ji and J() for the 
Neel and canted phases do not meet. We note, however, that simple (spline) extrapolations 
of the LSUBoo canted-phase magnetization curves generally give corresponding estimates 
for J2^ at which they meet their Neel-phase counterparts, which are in excellent agreement 
with those calculated as in Fig. 3. Such extrapolations also give clear evidence that the 
order parameter curves for the two phases meet at a value J2 = J2^ at which M is nonzero 
(and hence the transition is of first-order type) for all values of J[, and furthermore that the 
curves have no discontinuity in slope (or only a very small one) at J2 = ^2^. 

We remark that the one region where the extrapolation procedure for M becomes slightly 
problematic is for values of J[ < 0.1 near the (generally anisotropic) kagome point, J[ = 0. 
The reason is clear from Fig. 6(b) since the Neel curves drop to zero from a nonzero value with 
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a slope that approaches infinity as J( ^0. Thus, the Neel curves for the J[ = case extend 
(artificially) to the end-point J^"^^^^ ^ 0.592, whereas we know the value jlsuboo ^ q^^-^ ^^ 
this case. Clearly, the nature of the transition becomes rather singular in the limit J[ — » 0, 
as Fig. 6(b) clearly shows. 

We note finally that the CCM LSUBn solutions with n > 2 based on the canted state 
terminate at some upper critical value of J2 for all values of J[ (and Ji = 1). This provides 
preliminary evidence for another critical point at J2 = ^2^. For example, at the isotropic 
point J[ = Ji = 1 the termination points occur at values J2 ~ 74.5, 20.0, and 11.0 for 
LSUBn approximations with n = {4, 6, 8} respectively. To investigate this possible transition 
further we have performed a preliminary series of separate CCM LSUBn calculations based 
on the semi-striped state shown in Fig. 1(b). We find that this state is stable out to 
the J2 ^ 00 limit for all LSUBn approximations investigated {viz., with n < 8). The 
corresponding LSUBoo result for the semi-striped phase as J2 ^ 00 (for fixed J[ and Ji) 
is E/N ^ — O.2215J2, which may be compared with the exact result for this decoupled ID 
HAF chain hmit (of A^^ spins per chain) of E/{2N^) ^ -0.2216 J2. 

We note, however, that the LSUBoo results for E/N for the canted and semi-striped 
phases do not cross for J2 < 11.0, where the canted LSUB8 results terminate. Thus, it 
is not possible on this evidence alone to suggest that there might be a second first-order 
phase transition at J2 = J2^{Ji), with J2 = 1, between the canted and semi-striped phases. 
Nevertheless, the preliminary evidence is that the canted phase does not exist for values 
J2 > J2'{J[), with Ji = 1. 

V. DISCUSSIONS AND CONCLUSIONS 

In this paper we have used the CCM to study the infiucence of quantum fiuctuations on 
the zero-temperature gs properties and phase diagram of a frustrated spin-half HAF defined 
on a 2D square lattice with three sorts of antiferromagnetic Heisenberg bonds of strengths 
Ji, J[, and J2 arranged in the pattern shown in Fig. 1. The Ji and J[ bonds are between 
NN pairs on the square lattice, while the J2 bonds are between only those one-quarter of 
the NNN pairs shown. In the case when J[ = and J2 7^ Ji the model reduces to a spin- 
half HAF on a spatially anisotropic kagome lattice appropriate to the quasi-2D material 
volborthite. For the special case J[ = and J2 = Ji the model reduces to the spin-half 



15 



HAF on the (isotropic) kagome lattice. Similarly, for the special case J2 = and J[ = Ji 
the model reduces to the spin-half HAF on the (isotropic) square lattice. The model thus 
interpolates smoothly between these limiting cases. 

Classically the model has only two gs phases, namely an antiferromagnetic Neel phase 
and a ferrimagnetic canted phase shown in Fig. 1(a), and we have focussed attention in the 
present work on the effects of quantum fluctuations on these two classical phases. Consistent 
with the usual finding that quantum fluctuations favour collinear configurations of spins, we 
found that the phase transition point at J2^{J'i), with Ji = 1, between the Neel and canted 
phases satisfies the inequality J2^{Ji) > Jf = \{^ + ^^i) for all values J( 7^ 0, where J2 is 
the corresponding classical phase boundary, as shown in Fig. 4. Precisely at the anisotropic 
kagome-lattice point J( = 0, however, the classical and quantum critical values agree to the 
level of accuracy of our results; ^2^(0) = 0.51 ± 0.01, compared to Jf = 0.5. 

Our calculations provide strong evidence that the canted phase is not the stable gs phase 
for the model for values of J2 greater than some second critical value J2^{J[) with Ji = 1. 
Thus, unlike in the classical case where the canted phase is the stable gs phase for all 
values J2 > JfiJi) = 2^^ ~^ "^i^' ^^ ^^^ quantum spin-half case the canted phase seems 
to be the stable gs phase only for values ^2^(^i) < J2 < JTiJ'i)- ^"^ order to investigate 
the nature of the transition at ^2^(^i) with Ji = 1, we have also performed preliminary 
calculations using the semi-striped state of Fig. 1 (b) as model state in the CCM. Although 
we have convincing proof that such a semi-striped state is stable for the spin-half case under 
quantum fluctuations, for large values of J2 ^ 00 for all values of J[ and Ji = 1, its energy 
is always (very slightly) higher than that of the canted state in regions where solutions to 
the corresponding CCM LSUBn approximations both exist. While these preliminary results 
do not exclude a second first-order phase transition at J2 = ^2^(^i), with Ji = 1, from 
the canted phase to the semi-striped phase, it is also quite possible that the transition at 
J2 = J2^{Ji) is to an entirely different state. We hope to report further on the existence and 
nature of this second quantum phase transition in a future paper. 

As stated previously, our main aim here has been to discuss the entire phase boundary at 
J2 = Jo^iJi)-, with Ji = 1, of the model between the Neel and canted phases, for all values 
of the bond strength J[. Nevertheless, there is no doubt that the limiting case J[ = 0, 
corresponding to the spatially anisotropic kagome lattice is of huge interest in its own right, 
and we aim to discuss this case further in a separate future paper. The results for the onsite 
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magnetization M shown in Fig. 6 clearly indicate the special nature of the case J( = 0, as we 
discussed previously in Sec. IV. Figure 6(b) shows in particular that the order parameter M 
at J( = is a minimum for the isotropic kagome lattice (J2 = Ji = 1), and we have shown 
that our results for this isotropic case are compatible with the vanishing of the corresponding 
parameter M^ defined on the kagome lattice. Our results for the gs energy E/N for the 
isotropic kagome HAF also agree with the best available by other techniques (and see, e.g., 
Ref. [32]). 

The isotropic kagome HAF has been greatly studied in the past. The most direct re- 
sults from the exact diagonalization of finite lattices^^'^^ seem to give strong evidence for a 
spin-liquid gs phase. Such a conclusion is supported by block-spin approaches'^ '^^ and by 
various other studies. '''''^''^''*'^"^'* Nevertheless, conflicting results have been found by other 
authors^^"'^ who have proposed various valence-bond solid states as the gs phase of the 
isotropic kagome-lattice HAF. A detailed comparison of the exact spectrum of a 36-site fi- 
nite lattice sample of the isotropic kagome HAF against the excitation spectra allowed by 
the symmetries of various of the proposed valence-bond crystal states has, however, cast 
doubts on their validity. ^^ 

The classical ground states of the anisotropic kagome HAF are spin configurations that 
satisfy the condition that for each elementary triangular plaquette of the kagome lattice in 
Fig. 1 when the Bi sites and J[ bonds are removed (when J[ = 0), the energy is minimized. 
For J2 < I (with J[ = and Ji = 1) the classical ground state is collinear and unique, 
with the spins along the J2-bond chains aligned in one direction and the remaining spins on 
the kagome lattice aligned in the opposite direction. The total spin of this classical state is 
thus Stot = ^Nks where each spin has magnitude s. For the quantum case the Lieb-Mattis 
theorem^^ may also be used, for the limiting case J2 = only, to show that the exact gound 
state has the same value Stot = ^Nks of the total spin as its classical counterpart. 

For J2> \ (with J( = and Ji = 1) the classical ground state is coplanar with the canting 
angle = cos^^{jj- 7^ 0) shown in Fig. 1. The classical ensemble of degenerate coplanar 
states is now characterized by two variables for each triangular plaquette, namely the angle 
such that the middle spin of a given triangular plaquette forms angles (tt ± 0) with the 
other two spins of the same plaquette, and the two-valued chirality variable x = il that 
defines the direction (anticlockwise or clockwise) in which the spins turn as one transverses 
the plaquette in the positive (anticlockwise) direction. For a given value of J2 > ^ (with 
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J[ = and Ji = 1) the angle 7^ is given as above, and the different degenerate canted 
states arise from the various possible ways to assign positive or negative chiralities to the 
triangular plaquettes of the lattice. 

The HAF on the isotropic kagome lattice (with J[ = and J2 = Ji = 1) is especially 
interesting since for this case, with = |, the number Q of degenerate spin configurations 
grows exponentially with the number A'^^^ of spins, so that even at zero temperature the 
system has a nonzero value of the entropy per spin. By contrast, for the anisotropic case 
(with J[ = 0, J2 ^ Ji = 1), the degeneracy Q has been shown"^^ to grow exponentially with 



\/Nk [i.e., VL a exp(cyiVx)], so that the gs entropy per spin vanishes in the thermodynamic 
limit. Clearly, since in the limit J2 — ^ 1 the anisotropic model approaches the isotropic 
model, the anisotropic model must have an appropriately large number of low-lying excited 
states that become degenerate with the ground state in the isotropic limit, J2 ^ 1. 

The spin-half HAF on the spatially anisotropic kagome lattice has been studied by several 
authors recently using a variety of techniques. These have included large- A^ expansions of 
the Sp(A^) -symmetric generalization of the model, ^^ a block-spin perturbation approach to 
the trimerized kagome lattice, ^^ semiclassical calculations in the limit of large spin quantum 
number s,^^'^^ and field-theoretical techniques appropriate to quantum critical systems in 
one dimension (and which are hence appropriate here for the case J2 ^ Ji of weakly coupled 
chains). ^"^ The results of such calculations generally seem to indicate that the anisotropic 
kagome HAF (i.e., our model with J[ = 0) has a Neel-like gs phase, a canted coplanar gs 
phase and, in the limit of large anisotropy (J2 ^ Ji = 1), another gs phase that approaches 
the decoupled-chain phase as J2/J1 — * 00. The precise nature of this third phase is by no 
means settled, with the results of the various calculations not in complete agreement with 
one another. We hope to contribute our own more detailed CCM results to this debate in 
the two future papers outlined above. 
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